Polymer translocation through a nanopore under a pulling force 
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We investigate polymer translocation through a nanopore under a pulling force using Langevin 
dynamics simulations. We concentrate on the influence of the chain length N and the pulling force 
F on the translocation time r. The distribution of r is symmetric and narrow for strong F. We find 
that r ~ iV 2 and translocation velocity u ~ iV" 1 for both moderate and strong F. For infinitely 
wide pores, three regimes are observed for r as a function of F. With increasing F, r is independent 
of F for weak F, and then r ~ F~ 2+L ' for moderate F, where v is the Flory exponent, which 
finally crosses over tor~ F^ 1 for strong force. For narrow pores, even for moderate force r ~ F _1 . 
Finally, the waiting time, for monomer s and monomer s + 1 to exit the pore, has a maximum for 
s close to the end of the chain, in contrast to the case where polymer is driven by an external force 
within the pore. 

PACS numbers: 87.15.Aa, 87.15.He 



I. INTRODUCTION 

The transport of biopolymers through nanopores oc- 
curs in many biological processes, such as DNA and RNA 
translocation across nuclear pores, protein transport 
through membrane channels, and virus injection P, 
Due to various potential technological applications, such 
as rapid DNA sequencing [4j, gene therapy and con- 
trolled drug delivery la], recently a number of experi 
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and numerical 29. 
HI HE IH, El studies on polymer 
translocation have been carried out. 

In 1996, Kasianowicz et al. [6] demonstrated that 
single-stranded DNA and RNA molecules can be driven 
through the a-hcmolysin channel in a lipid bilayer mem- 
brane under an electric field, and the polymer length can 
be characterized. In addition, Li et al. [121 ] and Storm et 
al. [H[ showed that a solid-state nanopore could also be 
used for similar experiments. 



Inspired by the experiments [f .Il0l. 15] , a number of re 



cent theories [151, [17|, l20j, [27|, l28|, [29|, 130J, [31|, [33] have 
been developed for the dynamics of polymer transloca- 
tion. Particularly, the scaling of the translocation time r 
with the chain length N is an important measure of the 
underlying dynamics. Sung and Park [l7| and Muthuku- 
mar [201 ] considered equilibrium entropy of the polymer 
as a function of the position of the polymer through the 
nanopore. Standard Kramer analysis of diffusion through 
this entropic barrier yields a scaling prediction of r ~ N 2 
for the field free translocation. For the forced transloca- 
tion, a linear dependence of r on N was suggested, which 
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is in agreement with some experimental results ( 
for a-hemolysin channel. However, as Chuang et a, 
noted, the quadratic scaling behavior for the field free 
translocation cannot be correct for a self-avoiding poly- 
mer. The reason is that the translocation time is shorter 
than the equilibration time of a self-avoiding polymer, 
Tequii ~ N l+2u , where v is the Flory exponent [ill. HBl]. 
According to scaling theory, they showed that for large 
N, translocation time scales approximately in the same 
manner as equilibration time. For the forced transloca- 
tion, Kantor and Kardar [3(3] provided a lower bound 
for the translocation time that scales as N 1+l/ , by con- 
sidering the unimpeded motion of the polymer. Most 
recently, we investigated both free and forced translo- 
cation using both the two-dimensional fluctuatingbond 
model with single-segment Monte Carlo moves [32|, [33[ 
and Langevin dynamic simulations [34], [35[ . For the free 
translocation, we numerically verified t ~ N 1+2u by con- 
sidering a polymer which is initially placed in the middle 
of the pore [32|, [34| ■ For the forced translocation, we 
found a crossover scaling from r ~ N 2l/ for relati vely 
short polymers to t ~ N 1+v for longer chains [HI, |34| . 
In addition, we also found that this crossover scalin g re - 
mains unaffected for heteropolymer translocation [35j. 
The predicted short chain exponent 2v « 1.18 in 3D 
agrees reasonably well with the solid-state nanopore ex- 
periments of Storm et al. [TBI ]. 

Polymer translocation involves a large entropic bar- 
rier, and thus most polymer translocation phenomena 
require a driving force, such as an external electric field 
used in above mentioned experiments. However, one can 
also envisage the use of other forces, such as a pulling 
force. With the development of manipulation of single 
molecules, polymer motion can be controlled by optical 
tweezers [4j], [46|. This gives a motivation to study the 
translocation in which only the leading monomer expe- 
riences a pulling force. In addition, a new sequencing 
technique based on a combination of magnetic and op- 
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tical tweezers for controlling the DNA motion has been 
reported |47|. Therefore, it is of great importance to the- 
oretically investigate the polymer translocation under a 
pulling force. Kantor and Kardar [3(| have considered 
the scaling of r with N and with the pulling force F. 
They have also tested the scaling behavior in a MC sim- 
ulation study of the fluctuation bond model. However, 
as discussed below, this model is only valid for moderate 
pulling forces. For strong pulling forces the scaling of r 
with F is different and needs to be studied carefully. To 
this end, in this paper we investigate polymer translo- 
cation through a nanopore under a pulling force using 
Langevin dynamics simulations. 



II. MODEL AND METHOD 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains. The excluded volume 
effects and van der Waals interactions between all pairs 
of beads are modeled by a repulsive LJ potential as: 



U LJ (r) 



r > 
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(1) 



where a is the diameter of a bead, and e is the depth of 
the potential. Nearest neighbor beads on the chain is con- 
nected via a FENE spring with a potential U FENE{ r ) = 
— i/ci?Q ln(l — r 2 /i?g), where r is the separation between 
consecutive beads, k is the spring constant and i?o is the 
maximum allowed separation between connected beads. 

In the Langevin dynamics method, each bead is sub- 
jected to conservative, frictional, and random forces F; , 
Ff , and Ff , respectively, with [H] raf, = Ff +Pf +Ff , 
where m is the monomer's mass. Hydrodynamic drag is 
included through the frictional force, which for individual 
monomers is F^ = — £vj, where £ is the friction coeffi- 
cient, and Vj is the monomer's velocity. The Brownian 
motion of the monomer resulting from the random bom- 
bardment of solvent molecules is included through Ff 
and can be calculated using the fluctuation-dissipation 
theorem. The conservative force in the Langevin equa- 
tion consists of several terms Ff = — V(Ul,j + Ufene) + 
Fpuiiing- The pulling force is expressed as 



■ pulling 



Fx, 



(2) 



where F is the pulling force strength exerted on the first 
monomer and x is a unit vector in the direction perpen- 
dicular to the wall. 

In the present work, we consider a 2D geometry where 
the wall in the y direction is described as I columns of 
stationary particles within distance a from one another 
and they interact with the beads by the repulsive part 
of the Lennard-Jones potential. Wall particle positions 
do not change during the simulations. The pore is intro- 
duced in the wall by simply removing w beads from the 
wall. 



III. SCALING ARGUMENTS 
A. Stretching extension 

Here, we first briefly recall the main results of scaling 
analysis pioneered by Pincus [49[who considered a poly- 
mer under traction with two forces, F and —F, act on 
its end. The elongation L(F) of the chain maybe written 
as 



L(F)=i^(|), 



(3) 



where is a dimensionless scaling function, R = N v a, 
denotes the size of the unperturbed coil, and Q is the 
characteristic length of the problem, £ = ksT/F. 

For weak forces, such that F < ksT/N^a, the response 
is linear, i.e., 4>(x) ~ x, which leads to [44( 



L(F) ~ N 



2// 



k^T' 



(4) 



For moderate forces ksT/N^a < F < ksT/a, the chain 
breaks up into a one-dimensional string of blobs of size £. 
Then the elongation L(F) ~ JV. which leads to <p(x) ~ 



x 



ThuSj one obtains rug! 



kni 



(5) 



Finally for strong forces F > ksT/a, the chain is nearly 
fully extended with 



L(F) ~ Na. 



(6) 



In the following, we use similar arguments to analyze 
the extension of the tethered chain pulled with a constant 
force through a viscous medium. The geometrical imped- 
iment due to the finite width of the pore is neglected here. 
For clarity, we assume that the pulling force acts on the 
last monomer N. Without hydrodynamic interactions 
the force acting on segment n is given by 



(7) 



where Vi is the velocity of the ith segement. At steady 
state when inertia can be neglected compared with the 
frictional force, we assume v» = v and thus 



F n = n£,v. 



(8) 



Under this physical picture, we stress that the pulling 
force F equals Fn = N£v. Here, we encounter a situation 
in which the tension F n is segment-dependent and where 
the stretching of the chain is not uniform [5(J HH, [HJ . 

To this end, we generalize Eqs. (Op a nd ([5]) to this sit- 
uation, following Brochard-Wyart [50j], who considered 
the nonuniform deformation of tethered chains in uni- 
form solvent flow. Let ( n — ksT/F n be the n-dependent 
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size of the Pincus blobs and x n the position of the nth 
monomer with respect to the last monomer in the di- 
rection of the pulling force. For weak forces at n, i.e., 
F n < ksT/N^a, the local elongation at site n is 



dx r 



2v~l 

n a 



dn, 



(9) 



while for moderate forces, ksT/N^a < F n < ksT/a, 



k B i 



(10) 



Integrating Eqs. 
mation obeys 



and IjlOp over n, one finds the defor- 



L(F) ~ N 



2v- 



-l &<r 
k R T 



for weak forces and 



L(F) 



knl 



(11) 



(12) 



for moderate forces. From the scaling picture, the sizes 
£ n of the Pincus blobs obey ~ 1/F n ~ l/n ~ x~ v ' . 
This can be compared with the case of a tethered chain 
subjected to a uniform solvent flow [50j. The blob size 
decreases in the pulling force direction, resulting in a 
trumpctlike shape. For strong forces, the last part of 
the chain in fully stretched while its free end still shows 
Pincus-type behavior, corresponding to the regime called 
stem and flower [50] |. 

Note that Eqs. fTTj) and (fT2"|) can be obtained by sim- 
ply replacing F in Eqs. ([4]) and ([5]) by an effective force 
F e ff — N£v, which means Eqs.(|l|) and (O are still ap- 
plicable when one sets £ = £ e // = ksT/ F e ff . 



B. Translocation time 

To examine r as a function of N under the same con- 
stant pulling force F, we need to know L(F) as a func- 
tion of the pulling force F, not the drag force £v on 
monomer. To use Eqs. (jTTJ) and p2[) we need to use 
relation F = N£v, which are the same as Eqs.(H|) and (O 
although the microscopic pictures are different. 

Kantor and Kardar [30| have presented a scaling ar- 
gument for the unimpeded translocation, which serves 
as a lower bound for the true translocation through the 
nanopore. The polymer travels a distance L(i ? ))during 
the translocation process. The translocation velocity 
scales as F/N since the force is applied to one monomer 
only. Thus the translocation time should depend on 
N and F as r ~ 30]. For moderate forces i.e., 

ksT/N^a < F n < ksT/a, we have from the scaling of 
L{F) in Eq. flX 



L(F) 
v(F) 



jy2p-2+{ 



(13) 




FIG. 1: A typical configuration of a polymer of length N = 
300 pulled by a force of strength F = 5 during translocation 
process. 



This scaling relation for moderate force is the same 
as the one obtained earlier [3(J. We can now extend 
this approach to both weak and strong forces. For weak 
forces, the translocation time scales according to Eq.JTT]) 
as 



L{F) 
v(F) 



N 



l+2u 



(14) 



This scaling behavior is the same as that for translocation 
in the absence of forces [29|, [32|, H|| , in disagreement with 
recent claims [53| . For strong pulling forces, the polymer 
becomes fully stretched and thus obtains a qualitatively 
different spatial configuration. Such a configuration is 
shown in Fig. [T]for a polymer of length N — 300 during 
translocation. In this case, the translocation time scales 



as 



L(F) 
v(F) 



N 2 F~ 



(15) 



Eqs. (JT3J) and (fTS|l show that r ~ TV 2 for both moderate 
and strong forces. From Eqs. (fT4j) . (fl3|) and (fTS"]) . r as a 
function F have three regimes with increasing F. 



IV. RESULTS AND DISCUSSIONS 

In our simulations, the parameters are a = 1, fc^T — 
1.2e, and the time scale is given by t^j = (mcr 2 /e) 1 / 2 , 
which is in order of ps. The friction is set as £ = 
0.7m/t LJ . For the FENE potential, we use [H R = 2a, 
k = leja 2 . Typically, ksT/a — 4 pN for a chain with 
Kuhn length a = 1 urn at room temperature 295 K, 
and the time scale is about 11.28 ps for monomer mass 
to = 312 amu. The scale of the pulling force F is e/er, 
which is about 3.3 pN. The Langevin equation is inte- 
grated in time by a method described by Ermak and 
Buckholtz [55( in 2D. For the pore, we set w = 3er and 
/ = a unless otherwise stated. 

To create the initial configuration, the first monomer 
of the chain is placed in the entrance of the pore. The 
polymer is then let to relax to obtain an equilibrium con- 
figuration. In all of our simulations we did a number of 
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FIG. 2: The distribution of 2000 translocation times for a 
chain of length TV = 100 under the pulling force of strength 
F = 5. Here, the translocation times are normalized by their 
average value. 



runs with uncorrelated initial states. The translocation 
time is denned as the time interval between the entrance 
of the first segment into the pore and the exit of the 
last segment. The estimate for the translocation time 
was obtained by neglecting any failed translocation and 
then calculating the average duration of the successful 
translocations. Typically, we average our data over 1000 
independent runs. In this section, we only investigate 
the translocation under the consant pulling force. For 
the case with the constant pulling velocity, we will ex- 
amine it in detail in future. 



A. Translocation time distribution 

The distribution of translocation times for a polymer 
of length N — 100 pulled with a force F = 5 is presented 
in Fig. The histogram obeys Gaussian distribution. 
This distribution has a qualitatively different shape com- 
pared to that for the free translocation case, where the 
corresponding distribution is asymmetric, wider and has 
a long tail [29|, [Hj]. However, this distribution is quite 
similar to that for driven translocation under an electric 
field, in that it is narrow without a long tail and sym- 
metric [3(| IHj]. The stronger the pulling force, the nar- 
rower the distribution becomes. As a consequence of this 
distribution, the average translocation time r is well de- 
fined and scales in the same manner as the most probable 
translocation time. Of course, if a weak enough pulling 
force is used, we still can observe the long tail. 



B. Waiting time 

The dynamics of a single segment passing through the 
pore during translocation is an important issue. The non- 
equilibrium nature of translocation has a considerable 



effect on it. We have numerically calculated the waiting 
times for all monomers in a chain of length N. We de- 
fine the waiting time of monomer s as the average time 
between the events that monomer s and monomer s + 1 
exit the pore. In our previous work [33ll34| for transloca- 
tion under an electric field in the pore, we found that the 
waiting time depends strongly on the monomer positions 
in the chain. For short polymers, such as N — 100, the 
monomers in the middle of the polymer need the longest 
time to translocate and the distribution is close to sym- 
metric. However, for a polymer of length N — 300, it's 
approximately the 220 th monomer that needs the longest 
time to translocate on the other side of the pore. The 
waiting times for chains of length N — 100 and N = 300 
under pulling forces are presented in Fig . [3J As compared 
to the electric field driven case [33l l34j . the peaks of the 
waiting times are shifted towards the last monomers of 
the chain, independent of the force. This can be under- 
stood from the fact that when the chain is being pulled 
through the pore, its free energy increases due to loss of 
configurational entropy. For short chains, this leads to 
the chain motion to slow down almost monotonically un- 
til the chain entropies on both sides of the pore roughly 
balance each other. For long chains, the entropy of the 
pulled part eventually exceeds that of the remaining part 
of the chain and a maximum in waiting time appears close 
to the end of the chain. 



C. Translocation time as a function of the chain 
length 

As a reference point for comparison, we first check r 
as a function of N for an infinitely wide pore, as shown 
in Fig. EUa). We obtain in this case that r ~ /\ri-92±o.oi 
and r - /V 2 01±0 02 for F = 5 and 0.5 respectively. The 
value for F = 0.5 is in the moderate force regime while 
F = 5 corresponds to the strong force regime. For both 
these regimes, the theoretical prediction is t ~ N 2 , as in 
Eqs. (flU)) and (|15p . Our numerical results for these two 
pulling forces are in very good agreement with the scaling 
argument predictions. For a pore of finite width, the re- 
sults are shown in Fig. 0Jb). We get scaling exponents of 
1.87±0.01 and 1.91±0.02 for F = 5 and 0.5 respectively. 
These results are in excellent agreement with the Monte 
Carlo simulation results of Kantor and Kardar [3(| , and 
demonstrate that the scaling arguments for unimpeded 
translocation provides a useful estimate for the actual 
translocation through the finite size nanopore. 



D. Translocation time as a function of the pulling 
force 



Our theory predicts that there are three regimes in 
the dependence of the translocation time on the pulling 
force, as shown in Eqs. (TT4")) . (| 13[) and (|15p. To study 
this, we again consider first the unimpeded translocation 
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FIG. 3: Waiting times for a polymer of (a) length N = 100 FIG - 4: The translocation time as a function of polymer 
with F — 5, (b) length N = 300 with F = 5. len S th N for ( a ) an infinitely wide pore and (b) a pore of 

finite width. A constant pulling force of strength F — 0.5 
and 5 acts on the first monomer. 



through an infinitely wide pore. The numerical results in 
Fig. E^a) confirm the existence of the three regimes. The 
translocation time is independent of F for weak pulling 
forces, i.e., F < 0.3, which is indicated in Eq. (fFi| . With 
increasing pulling force, the translocation time scales 
with the force with an exponent of -0.67 for 0.3 < F < 2. 
This result is in good agreement with the theoretical pre- 
diction in Eq. {T®, where r - F" 2 +^ - in 2 d. 
However, for 2 < F < 10, the exponent is -0.84. This 
shows that L(F) ~ ( ^p^)" 1 is no longer valid since we 
are in the strong force regime of F > ksT/a. Instead, 
we expect L(F) to be almost independent of the force 
and correspondingly the translocation time should scale 
as Eq. (fl~5|) in the limit of a strong force. 

For a pore of finite width, the translocation time as a 
function of pulling force is presented in Fig. [5jb). For 
0.25 < F < 10 we find r - ^-o.94±o.oi^ Thege results 

show that under the restriction of the wall for a pore of 
finite width, it is much easier for the chain to become 
fully stretched and hence the strong force limit scaling 
behavior holds r ~ F^ 1 through the entire range of forces 
applied. 

Finally, it is important to note that both in the case 
of external field (voltage applied across the pore) [12, [HI 



and pulling force driving the translocation process, there 
exists a fundamental difference between the Monte Carlo 
results for the lattice fluctuating bond model and the 
continum model considered here in the strong driving 
force limit. In the Monte Carlo study, the microscopic 
transition rate saturates very quickly when the external 
driving force increases, leading to a saturation of the ve- 
locity and the translocation time [13, H3|- This aspect 
of the fluctuating bond model is unrealistic and does not 
correspond to the true dynamics of the system. The con- 
tinum model does not suffer from this artifact. As seen 
in Fig. [51 in the present model the translocation time r 
scales asr~ F^ 1 up to the maximum force value studied 
and shows no sign of saturation. Our previous studies of 
the field driven translocation with both Monte Carlo and 
Langevin dynamics show that while the scaling behavior 
agree in most regimes, the same difference occurs in the 
strong force limit. 
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CONCLUSIONS 
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In this work, we have investigated the dynamics 
of polymer translocation through a nanopore under a 
pulling force using 2D Langevin dynamics simulations. 
We have focused on the influence of the length of the 
chain N and the pulling force F on the translocation 
time r. The distribution of r is symmetric and narrow 
for strong F. We find that t ~ N 2 and translocation 
velocity v <~ TV -1 for both moderate and strong F. For 
infinitely wide pores, three regimes are observed for r as 
a function of F. With increasing F, r is independent of 
F for weak F, and then r ~ F~ 2+v for moderate F, 
where v is the Flory exponent, which finally crosses over 
to r ~ F^ 1 for strong force. For narrow pores, even for 
moderate force r ~ F^ 1 . Finally, the waiting time, for 
monomer s and monomer s + 1 to exit the pore, has a 
maximum for s close to the end of the chain, in contrast 
to the case where polymer is driven by an external force 
within the pore. 
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